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_. Abstract 

ON 

In 2011 Radyushkin outlined the practical implementation of the One-Component 

, , Double Distribution formalism in realistic Generalized Parton Distribution models. We 

rG compare the One-Component Double Distribution framework to the standard one and 

compute Deeply Virtual Compton Scattering observables for both. In particular the new 
implementation is more flexible, offering a greater range of variation of the real and imag- 
inary parts of the associated Compton Form Factor while still allowing to recover results 
similar to the classical approach. Moreover the polynomiality property is satisfied up to 
,__! the highest order. Although the comparison to experimental data may be improved, the 

One-Component Double Distribution modeling is thus an attractive alternative. 
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Introduction 

m 

Generalized Parton Distributions (GPDs) were introduced by Miiller et al. [lj, Ji [2J and 
Radyushkin |3J. They encode a wealth of information about the structure of the nucleon 
including 3D imaging of its partonic content and access to the quark orbital angular momentum. 
GPDs have been the object of an intense theoretical and experimental activity ever since (see 
the reviews Ref . [H El El [8j [9] and references therein) . 

GPDs can be accessed experimentally by studying the processes of leptoproduction of a 
photon: Deeply Virtual Compton Scattering (DVCS) (GPDs can also be accessed through its 
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timelike counterpart: Timelike Compton Scattering [10J), or a meson: Deeply Virtual Meson 
Production (DVMP). The DVCS amplitude is expressed in terms of Compton Form Factors 
(CFFs), which are convolutions of GPDs with known kernels. CFFs were extracted from DVCS 
in Ref. P^l H21 H31 QH HSl HS1 HZl US] while GPDs were constrained from DVCS in Ref. [32] or 
from DVMP in Ref. [HUSH [22]. Note that this last set of GPDs, tuned for DVMP analysis, has 
been compared recently to almost all existing DVCS measurements [23J . Although these results 
are promising, our knowledge of GPDs is far from complete and will certainly be challenged 
by forthcoming measurements at COMPASS, by the high precision data expected after the 
Jefferson Lab 12 GeV upgrade and at an Electron Ion Collider. Gaining an experimental 
knowledge of GPDs is more involved than extracting Parton Distribution Functions (PDFs) 
from measurements: there are more GPDs than PDFs, and they depend on three variables 
instead of one. GPDs are subject to several theoretical constraints, and are related to PDFs 
and quark contributions to nucleon Form Factors (FFs). There is no known parameterization 
of GPDs relying only on QCD first principles. 

Part of the theoretical constraints on GPDs are fulfilled by modeling Double Distributions 
(DDs) PJ I2U [25], which are related to GPDs by a Radon transform [26]. DD modeling is the 
most popular way to build simple and realistic GPD models, and has been widely used so far, 
for example in the popular VGG (Vanderhaeghen, Guichon and Guidal) 0, E3 EH [291 EH] or 
GK (Goloskokov and Kroll) [201 I2T1 [22] models. However these models are not in agreement 
with all existing DVCS data, as one can see from fit results for VGG [LLJ 1X2} [T3j fH] or in 
the systematic application of the GK model to DVCS data in Ref. |23j . To reach a better 
comparison to the data one may either turn to more sophisticated GPD parameterizations, for 
example as advocated in Ref. |31| . or try different implementations of DD modeling. The latter 
direction recently benefited from renewed interest when Radyushkin exhibited a GPD model 
[32] built on the so-called One- Component DDs representation [7j. This representation was 
first outlined in Ref. [33] but was not used for phenomenology so far: its implementation was 
an open problem since it leads to a more divergent behavior of the GPDs at small longitudinal 
momentum fractions than the one described in Ref. |2"^[ 123] . 

In this paper we implement the One-Component DD formalism along the lines described in 
Ref. [32] for PDFs with an integrable singularity at the origin, and confront it to DVCS experi- 
mental data. We will compare our model to Jefferson Lab Hall A [31] and CLAS measurements 
[35] which are very accurate and concern the valence sector. Therefore we only implemented 
the One-Component DD model for the singlet contribution of valence quarks. Since the par- 
tonic interpretation of DVCS observables relies on factorization theorems, we will apply the 
further restriction \t\/Q 2 < 0.1 where t is the momentum transfer, and Q 2 the incoming photon 
virtuality [3j [361 ETJ [38]. The work of Ref. [32] was done in the case of a spinless target although 
it was illustrated by a nucleon valence PDF-like toy model. The analog of the GPD H spin0 in 
the spinless case is an admixture of the GPDs H spin i/ 2 and -E'spini/2- As discussed in Ref. [23], 
beam polarized DVCS observables involving unpolarized targets are mostly independent of the 
GPD E. This allows us to use the results of Ref. [32] without any further adaptation to the 



spin- 1/2 case in this first study. 

In the first section we remind the basics of GPD modeling based on DDs. In the second 
section we modify accordingly the valence part of the DD model used in Ref. [23j and in the 
third section we discuss the phenomenological consequences of this alternative DD Ansatz. In 
the fourth section we discuss the model for nucleon GPDs relying on both 1CDD and 2CDD 
formalisms recently advocated in Ref. [39J. 

1 GPD models in the Factorized Double Distribution ap- 
proach 

For any four- vector v we note: 

v ± = —^(v° ±v 3 ) and v = {v + ,vj_, v~). (1) 

v2 

(uv) denotes the scalar product of two four-vectors u and v. 

1.1 One- Component DD and Two-Component DD formalisms 

For the sake of clarity, we discuss in this section the case of GPDs and DDs of spinless hadrons. 
The GPD H q , q denoting the quark flavor, is introduced through the matrix element of 
Eq. @: 



1 f dz ixP +~- I ~ A 



H q (x,^t) = - / e lxF z ~ { P , 

v '*" ; 2 / 2tt \ 2 



Z\ , I z 

78 



z i =0 



where ^ = — A + /(2P + ) is the skewness. 

The DDs F q and G q of the two-component DD (2CDD) formalism associated to the quark 
flavor q are defined by the following matrix element: 
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H) 7 ^(i)| P -f) 20 = ^l^ae-^^^F q ^,a,t) 

-A M I d/3dae- i ^ Pz)+ia ^ 1 G q {f3,a,t) 
Jn 

+ higher twist terms, (3) 

where q and q denote the quark fields separated by the light-like distance z, P ± A/2 the 
momenta of the incoming (+) and outgoing (-) hadrons and t = A 2 is the usual Mandelstam 
variable. When not necessary, the t-dependence will not be explicitly written. Q is the rhombus 
defined by: 

\a\ + \P\<l. (4) 



This yields the following relation between GPDs and DDs: 

H q (x, £,t)= [ d/3da 6{x-/3- af) (F 9 (/3, a) + £G 9 (/3, a, £)) . (5) 

Ju 



Assuming the vanishing of DDs on the boundary of the rhombuqjfi, one can see that gauge 
transformingn [2GJ 0U] the DDs F and G by means of an arbitrary function a q : 

F q (P,a) -> F*(J3,a) + ?£(J3,a), (6) 

<?«(£,«) -> G^,a)-^(Aa), (7) 

does not change the matrix element on the left-hand side of Eq. rt3j), neither does it modify 
the resulting GPD in Eq. ([5|. From time reversal invariance G 9 and a 9 are a-odd while F q is 
ct-even. 

A representation relying on one unique DD f q ((3, a) was proposed in Ref. [33J and is named 
One-Component DD (1CDD) in Ref. [7J. The non-trivial (5 and a dependence of the F and G 
type DDs is expressed in terms of / by the following: 

*icdd(M = /W,a), (8) 

GicddGM) = <xf q {P,a)- (9) 

This DD representation is less known than the 2CDD but is equivalent as it can be obtained 
from the 2CDD representation by a gauge transform (see Ref. [10]). In the 1CDD representation 
Eq. pi can thus be written: 
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HH g (i) | P "f ) 2 = l^ae-^ p ^^(2P^-A,a)f q ^a,t) 

+ higher twist terms. (10) 



With our forthcoming application to DVCS data in mind, let us focus on singlet GPDs and 
DDs, defined by: 

H q+ (x,0 = H q (x,0-H q (-x,0, (11) 

F q+ (f3,a) = F q ((3,a)-F q (-(3,a), (12) 

G q+ (/3,a) = G q (/3,a) + G q (-/3,a). (13) 



4 As stated in Ref. [40 DDs need only vanish at the corners (/3 = ±l,a = 0) and ()3 = 0, a = ±1) of the 
support 51. This fact introduces an extra boundary condition which does not modify our discussion and thus is 
omitted. 

5 The terminology of gauge tranformation has been chosen in Ref. |5B] for its formal likeness with the gauge 
transformation of the vector potential of the static two-dimensional magnetic field. 



The /3-odd <Td function: 
(T q D (f3,a) = -- 



dp'G q+ (P',a)- / df3'G q+ ((3', a) - sgn(/3)D q (a) 

l+[or| i/3 



(14) 



reduces the G-type singlet DD G^ D+D (/3, a) to a function 5([3)D(a) with a trivial dependence 
on the variable (3: the so-called Polyakov - Weiss D-term |H] while all the j3 dependence goes 
to the F-type singlet DD F^ D+D ((3,a). This representation is coined "DD+D" in Ref. |32| and 
corresponds to the Polyakov - Weiss gauge. 

A popular way to model the DD F q is Radyushkin's Factorized Double Distribution Ansatz 
(FDDA) |l2] which makes contact with the forward limit q(x) of the GPD H q (x) when t = 0, 
or with the (Fourier transformed) unintegrated forward limit q(x, t) when t ^ 0: 

F q (/3,a,t) = 7T N (/3,a)q(P,t), (15) 

where the profile function n N reads: 

MP > a) = ^Y{ N+ i) (i-i/3ir +1 • (16) 



r denotes the Euler Gamma function. The coefficient T(N + 3/2)/(^ ] /^TT(N + 1)) guarantees 
the normalization of the profile function: 

da7TAr(/3,a) = 1 . (17) 

1+1/31 

From a matter of principles, the FDDA can be applied to a F-type DD in the 1CDD and 2CDD 
formalisms as well. However the FDDA breaks the equivalence between the two representations 
as we can see by implementing the gauge transformation between the DD+D and the 1CDD 
representation: from Eq. (pi) we see that: 

do I da = ^ CDD (/3,a)-^ D+D (/3,a) 
= n N (/3,a)q((3) -n N ((3,a)q(/3) 
= 0. (18) 



and a does not depend on a; Being constant and a-odd, a must be zero. But o~d in Eq. (|14j) is 
manifestly non-trivial. Choosing to apply the FDDA to one representation or another is thus 
a phenomenological choice and can a priori be decided by comparison to experimental data. 

One of the advantages of the 1CDD parameterization involves the polynomiality of Mellin 
moments: 

/ dx x n H q (x,0 = X>0 1 < +1 , + mod(n,2)(20" +1 ^ +1 , (19) 

J —1 t—n 



where the scalar coefficients A q n+1 i and C^ +1 depend on t. Indeed, in the 2CDD (DD+D) 
representation, the term of highest degree is generated by the G-type DD (D-term), i.e. the 
F-type DD alone does not fulfill the polynomiality property. On the contrary, in the 1CDD 
representation, all powers of £ are correctly generated by one DD only. 

Note however that the FDDA applied to the 1CDD formalism leads to more singular GPDs. 
Indeed, in that case the factorized Ansatz yields: 

fo(P,a) = q M M ^ a y (20) 

The relation in Eq. d5J) between the GPD H q and the DDs in the 1CDD framework (see Eq. (J8J), 



Eq. Q and Eq. gjf) reads: 

H\x,0 = [ dj3da5(x-j3-at)(F? CBB (j3,a)+£Gl CDD (P,a)) 

Jn 

= f df3da5(x-f3-a£)(/3 + a£)f q (l3,a) 
Jn 

= x I d/3daS{x-/3-a^)^y-Tt N {^,a). (21) 

Jn P 

Given the typical behavior of nucleon valence PDFs q va \((3) oc /3~ 05 for small /3, the divergence 



at /3 = in Eq. (21 ) is problematic. Looking at sea quarks, the divergence becomes even worse, 
as the behavior of the associated PDFs is given by q se n((3) oc /3~( 1+<5 ) with < 5 < 1, when (3 
goes to 0. 

To unambiguously define the valence and sea contributions to PDFs and GPDs, we follow 
the conventions of Ref. |43| and references therein: 



• The valence and sea contributions g va i and g sea to the PDF q are defined on [—1, +1] by: 

<Zval(/3) = 0O%val|[O,l](/?), (22) 

g sea (/3) = sgn(/3)g sea | [0i i](|/3|), (23) 
where <Z|[o,i] denotes the restriction of the PDF q to the interval [0, 1]. 

• The valence and sea contributions _F^ al and F^ ea to the DD F q are: 

*&(&«) = (i^(/3, «) + F«(-/3, «))#(/?) (24) 

if ea (/3,a) = (F\^a)6{P)-F\-^a))6{-P) (25) 



To take care of this decomposition the classical FDDA (15) writes: 

FUP,ot) = ic N ^(P,a)<h*(P), (26) 

F^ifra) = 7r^ ea (/3,a)g sea (/3), (27) 

where the profile function parameters A^ va i and A sea can be distinct numbers. 



1.2 Taming the divergent structure of the One- Component DD for- 
malism 

In this section we describe Radyushkin's solution 



to control the potential divergence 
of Eq. (21) and prepare the ground for an implementation into a realistic GPD model. More 



precisely, Radyushkin made manifest that the 1CDD structure of the GPD model of Szczepaniak 
et al. derived in Ref. |45| . and explained how the aforementioned divergent behavior at small 
parton longitudinal momentum fraction can be controlled through a dispersion relation. 



k x = k- A/2 



k 2 = k + A/2 




P-A/2 P + A/2 

Figure 1: Quark hadron scattering amplitude. 



The considered model is essentially the computation of a triangle diagram. It relies on three 
main assumptions (see Fig. II] for notations): 

(i) A parton - hadron scattering amplitude T(P, A, k) in replacement of the spectator quark 
propagator. 

(ii) A once-subtracted dispersion relation for the quark-hadron scattering amplitude T(P, A, k): 



T(P,A,k)=T (t)+ / dap{a 



+oo 



a-(P-k) 2 a 



(28) 



where the independent subtraction constant T is unknown and p is a spectral function 
that generates a Regge-behavior for PDFs. 



(iii) A modification of the quark propagators: 

1 1 



— > 



(m 2 - jfe?) (m 2 - kj) N 



for i = 1,2. 



(29) 



The spectral function p can be traded for the forward limit q va \(x) of the GPD H^ al (x, 



KJ^O 



x 



d/3 



+1 -" qUP) ,* , 
da — - — vtat(p, a) 

1+13 



P 



S(x — (3 — a£) 



(l-/3) : 



;S X 



1=0 



Having remarked that ttn ( (3, a(l — /3)j = vrjv(0, a)/(l — f3), a change of variables yields: 



(30) 



x 



dp 



+1-/3 



da5(x — (3 — a£) 



1+/3 



Qval{P) , . . . , . f Qyal{j) 

— - — ir N (/3,a) -6{l3)Tc N (0,a) / d7 



This last equation clearly displays the 1CDD structure of this model, with: 



o 7(l-7) 2 
(31) 



(32) 



The singular behavior of the forward function q ya i(P)/P is regularized by the second term 



in the brackets of Eq. (30) as will be shown explicitly below in Eq. (41) and Eq. (42) when 



implementing the model. This term stems from the 1/er term of the dispersion relation Eq. (28). 
At last, the subtraction constant T Q generates an additional term D^\ 



mt) x 



D °\£'V 2N2 2N £ 



x 

e 



N 



which can be recast in the 1CDD framework by introducing: 



\x\ <£), 



D q Ja) 



(33) 



/ 9 (/3,a)- / dpda6(x-0-ag)8(B)-± 

In " 



(34) 



where the variable t has been omitted for clarity. 

2 Implementation 

As stated in the introduction, the VGG and GK models are probably the DD models that 
are most often used for phenomenological applications. Both are expressed in the 2CDD for- 
malism, in the specific DD+D representation. Moreover the D-term is set to in the GK 
implementation. 

The GK model was built in order to interpret the experimental results on Deeply Virtual 
Meson Production (DVMP), and as such, it was designed to work at small to intermediate Xb 
values. It was recently confronted to DVCS measurements in Ref. [23J. Without any tuning 



6 See Eq. (16) of Ref. [44] ; there is a typing mistake in Eq. (36) of Ref. [32 modifying the expression of the 
regularizing term S[x — a/(l — /3)£]. 
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of parameters, the model reaches a good agreement to the data at small xb i-e. in the sea 
quark region. Not surprisingly the comparison is less satisfactory in the valence region, leaving 
open the question of the extension of the model's validity range. A natural question is: Can 
the agreement be improved in the valence region by switching from the 2CDD to the 1CDD 
formalism? 

2.1 Valence part of the GPD H in the Goloskokov - Kroll model 

The valence GPD i^ al for the quark flavor q is described in the 2CDD formalism, in the DD+D 
representation, assuming a vanishing D-term. In the following this specific representation will 
be simply referred to as "DD". A factorized Ansatz is used: 

H^{x, £, t, fi 2 ) = f d(3da tt n (P, a)9((3)q val ((3, t, fi 2 )6(x - - a£). (35) 

Jo, 

The ^-dependent PDF q va ,\((3,t, /i 2 ) is parameterized as: 

2 

g va i(/U/i 2 ) = /3- Q V 5 (l-/3) 2n+1 E c ^ § - ( 36 ) 

3=0 



The coefficients 5 and Cj have been determined in a fit to the CTEQ6m PDFs [46J with n = 1. 
This specific choice of the PDF parameterization allows the analytic computation of the GPD 
Hl &l . In practice, we computed numerically the GPD H^ and checked that the analytic 
calculation gave the same result. The coefficient a' has been chosen in order to approximate 
the small ^-dependence of the quark contribution Fl to the proton form factor F\. The values 
of these coefficients for u and d quarks are recalled in Tab. [T] 

In Ref. |20l I2T1 1221 123] the profile function exponent iV is fixed to 1, but in the following we 
will allow it to take real values between 1 and +oo. Such a change of N does not modify the 
description of the form factor F\\ 

Ff{t) = f dx /d/3davr iV (/3,a)^(/3)g val (/3,t,/i 2 )(5(x-/3-aO 
J-i Jq 

1 /■+!- 



/ d/3g val (/3,t,/i 2 ) / da7rjv(/3,a) (37) 

JO J -1+0 

L 

d/3g val (/3,t,/i 2 ). (38) 



In this derivation we used the normalization (17) between Eq. (37) and Eq. (38). 
Using the following notation: 

h = ^4 and ft = ^£ (39) 





Mval 


dva.1 


5 


0.48 


0.48 


Co 


1.52 + 0.248 L 


0.76 + 0.248 L 


C\ 


2.88 - 0.940 L 


3.11 - 1.36 L 


c 2 


-0.095 L 


-3.99 + 1.15 L 


a'(GeV- 2 ) 


0.9 


0.9 


n 


1. 


1. 



Table 1: parameterization of the valence quarks PDFs in Eq. (36) with L 
Ql = A GeV 2 . 



log(QVQo) and 



Eq. (35) can be further simplified: 

1 






(40) 



The dependence of the PDF on the factorization scale fi 2 is approximated through the in- 
dependence of the PDF coefficients exhibited in Tab. [TJ In the GK model, this is the only 
dependence of the GPDs on the factorization scale. The factorization scale is chosen to be 
equal to the photon virtuality : fi 2 = Q 2 . 



2.2 Modification of the valence sector 

In Ref. [32J Radyushkin used a nucleon-inspired PDF toy model q(/3) = (1 — (3) 3 /y/]3 in the 
1CDD framework. The resulting GPD is peaked at x ~ £, the height of this peak being a 
decreasing function of N. This means that the imaginary part of the Compton Form Factor l-i q 
computed at Leading Order (LO) in QCD perturbation theory (i.e. Im'H LO (C,) = it H q+ '(£ , £)) 
may depend markedly on the profile function exponent N. This is of direct phenomenological 
relevance since measurements in the valence region with a polarized beam directly depends of 
this quantity [4711481149]. 
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Following Ref. |32j, the analog in the 1CDD framework of the 2CDD expression Eq. (40) is: 
H*4(\x\<te,t,lt) 



X [?* q val (f3,t,fi 2 ) 



*„[0,^)-ic N UZ(i-p) 



X [ h gyajQM^ ( R X 



i - 



i 



1-/3 



d/3 






Hl l (x>Z,£,t,fi 2 



x f^ AR qU^ty) f R x-/3 



(41) 
(42) 



Starting from Eq. (30) the F\ form factor sum rule writes: 

■+i 



Ff{t) 



dxx [ d(3dan N (P,a)6(P)^^^ 
i Jn P 

5(x — f3 — a£) — — 7^5 I x — 



+i 



o 



■+1-/3 

d/?gvai(/3, t, //) / da7rAf(/3,a) 

1+/3 



Jo 

+1 



/? 



1- 



1 



(1-/3) J 



+1-/3 



daa7TAr(/3, a) 



1+0 



d^g va i(/5,t,/i 



(43) 
(44) 







In this derivation, we used the normalization (17) and the fact that n^ is an even function 



of a between Eq. (43) and Eq. (44). The lessons of this computation are twofold: firstly, the 



result Eq. (44) is the same as Eq. (38). This means that we can change the DD formalism from 



1CDD to 2CDD without spoiling the t-dependence tuned to fulfill the form factor sum rule. 



Secondly, the final expression of the sum rule in Eq. (44) does not depend on the exponent N 



of the profile function. Thus we can also change N without altering this sum rule. Then to 
observe the phenomenological consequences of the choice of the DD formalism, we just have to 



change Eq. (40) to Eq. (41) or Eq. (42). 



When iV — > +oo the profile function iTN(/3,a) — > 5(a). Taking this limit in the classical 

9(x)q va i(x,t,fi 2 ). 



2CDD formalism, Eq. (35) yields: 

HUx,H,t^ 2 



Noo 



(45) 



Following the terminology of Ref. [50], we will quote the resulting asymptotic GPD model 



Forward Parton Density (FPD). The same exercise applied to the 1CDD model Eq. (30) also 
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gives: 



#4iO>£>^ 2 



Noo 



6{x 



g va i(a;,t,/i 2 



x 



(46) 



where we used x5(x) = 0. When x 7^ both 1CDD and 2CDD models have the same limit at 
large N. The resulting GPDs do not depend on the skewness £. 

Consequently, one may ask how the skewness ratio i.e. H(x, x)/H(x, 0) at small x is mod- 
ified. Using a Regge-behaved PDF toy model q(x) oc x~ a , one can compute this ratio analyti- 
cally in the small x limit. It has been done for example in Ref. [17] for the 2CDD formalism. 
We computed this ratio in the case of 1CDD formalism for < a < 1 because we used the 
once-subtracted dispersion relation Eq. (|28J). Both results are different: 



H q (x,x) 



Hi(x,0) 
W(x,x) 



Hi(x,0) 



2CDD 



1CDD 



y _ a r(2iv + 2) r(iv + i 
' r(jv + i) 

t . a T(2N + 2) 



a) 



T(2N + 2-a) 
T(N + l-a) 



2N + l-a 



T(N + 1) T(2N + 2-a) 2(N - a) 



H«(x,x) 



Hi(x,0) 



2N + 1 



a 



2CDD 



2(N-a) 



(47) 
(48) 
(49) 



The evaluations of the skewness ratio go to 1 as iV grows to infinity. On Fig. [2] we compare the 
dependence on iV of those two ratios at an arbitrary chosen value of a g]0, 1[. When both a 
and N go to 1, the skewness ratio is divergent in the 1CDD case. 




Figure 2: Skewness ratio at the arbitrary value a = 0.5 as a function of the profile function 
parameter N for the 1CDD (solid red line) and 2CDD (dashed blue line) parameterizations. 
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3 Comparison to experimental data 

As stated before we will apply this 1CDD GPD model to DVCS measurements in the valence 
region. We work with a Leading-Order (LO) definition of the CFF 7i q : 

i / 1 1 



H g (S,t) = / dxH%x,t,t)[t r-- -), (50) 

ReU q (£,t) = pj dxH«(x,£,t)(^-j^A, (51) 

ImW(1-,t) = 7r(#*(z,£,i)-iF(-x, £,*))■ (52) 

where "P denotes Cauchy's principal value prescription. We denote "H the average of "H^s 
weighted by the square of quark electric charges. 

It is commonly believed that Next-to-Leading Order (NLO) corrections have a small impact 
in the valence region, hence justifying the LO approximation. However complete next-to-leading 
order expressions of CFFs are available (37J EH |52j [53j EPJ Ell ESI EH EI] and recent estimates 
[58J challenge the aforementioned common view: quark and gluon NLO contributions may not 
be negligible even in the valence region at moderate energy. These results triggered an ongo- 
ing theoretical effort on the soft-collinear resummation in DVCS |59| 150] . New developments 
in this direction are expected in the near future. However our concern here is not a detailed 
phenomenological study of DVCS in the valence region, but rather the test of a new param- 
eterization of the GPD H . Therefore we will stay close to the prescriptions used in recent 



evaluations of DVCS observables [23] and use the LO expressions Eq. (51) and Eq. (52) of the 
CFF H. 

As stated in the introduction, we will apply the 1CDD formalism to the computation of 
JLab Hall A and CLAS polarized beam observables, namely beam helicity dependent and 
independent cross sections [34J and beam-spin asymmetries [35] . All details about the evaluation 
of observables related to the ep — > ep'-f channel are given in Ref. |23j . In particular we use the 
Trento convention [61j for the definition of the angle between the hadronic and leptonic planes. 
The variable xb is classically defined as: 

q 2 

x B = —z , (53) 

2p- q 

where p and q are the 4-momenta of the target nucleon and the exchanged virtual photon in 
the Born approximation. We also take the following asymptotic formula for £: 

f = -^-- (54) 

^ 2-x B K ' 

When comparing our model with the data, one should keep in mind every limitation it 
contains, including the leading order in a s approximation, the modification of the valence sector 
alone, the neglect of higher twist corrections or the treatment of QCD evolution equations. 
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Figure 3: Handbag diagrams for the DVCS process. 

3.1 Tuning of the profile function 

In the 2CDD formalism it is common to set the exponent iV = 1 in the profile function n^ 
because it corresponds to the asymptotic shape of a non-singlet quark distribution amplitude; 
this is the choice made in the GK and VGG models. However here we study DVCS data with 
typical Q 2 values between 1 and 4 GeV 2 , and xb between 0.1 and 0.5. It is not clear at all 
that the asymptotic shape is reached. It is thus interesting to compute the CFF % and DVCS 
observables for different values of the profile function exponent N. 



In a first step, we fit our GPD parameterization Eq. (41) and Eq. (42) to Hall A beam 
helicity-dependent and independent cross section data jM]. Although measured on a restricted 
kinematic domain, these data are highly accurate, and as emphasized for example in Ref. [23J, 
the helicity-dependent cross sections are mostly sensitive to the singlet GPD H + (£,!;, t,Q 2 ). 
This offers the possibility to set the width of the profile function through the parameter N, and 



independently to quantify the impact of an additional subtraction term such as Dq in Eq. (33). 

First we compute the CFF % obtained in the 1CDD and 2CDD pictures at the given 
kinematic configuration xb = 0.36, Q 2 = 2.3 GeV 2 and t = -0.23 GeV 2 , and compare it to the 
CFF obtained in the N — > +oo limit (FPD model). We were able to compute CFFs with iV 
ranging from 1 to 50 and managed to get the correct asymptotic behavior, which is a test of 
our numerical integration routines. On Fig. [4] we observe that the 2CDD model with N = 1 
is almost equal to its asymptotic limit. This fact has been known for some time [5]. However, 
we see a marked difference with the 1CDD implementation for small values of N. Therefore, 
we expect a clear difference when computing actual observables. We also expect that for large 
enough N, the 1CDD parameterization will produce similar results to the 2CDD formalism 
since both have the same asymptotic limit. 

Hall A helicity-dependent and independent cross sections, restricted on kinematics such that 
\t\/Q 2 < 0.1 , are best described by 1CDD choosing N ~ 1.86 (see Tab. ^1). The comparison 



to part of these data is shown on Fig. p In this comparison, all GPDs (H, E, H and E) are 
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Figure 4: CFF T-L as a function of the exponent N of the profile function tcn at x# = 0.36, Q 2 
= 2.3 GeV 2 and t = -0.23 GeV 2 for the 1CDD (solid red line), DD (dashed blue line) and FPD 
(dash-dotted green line) parameterizations. 



taken into account in a leading-twist leading-order evaluation along the lines of Ref. [23j. In 
particular the valence and sea parts of the GPD H are computed, but only the valence part 
has been changed from the 2CDD to the 1CDD formalism. We observe a discrepancy between 
the 1CDD prescription and beam helicity-dependent cross sections. This is due to the fact 
that these data are consistent with the FPD model which is formally accessed in the 1CDD 
formalism when N — > oo. From Fig. [4] we see that the 1CDD CFF H at N ~ 5 is already a 
good approximation of the FPD CFF H. Moreover, the smaller iV is, the larger IrriH becomes 
in the 1CDD picture as displayed above in Fig. |4} On the contrary, this 1-parameter fit leads 
to a much better agreement with helicity-independent cross sections for (f) close to 0° or 360° 
and 6 close to 180° as well. 
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Parameterization 


rVd.o.f 


1CCD (N ~ 1.86) 

DD 

FPD 


4.0 
5.9 

7.2 



Table 2: % 2 per degrees of freedom for the comparison of the different parameterizations with the 
subset of Hall A beam helicity-dependent and independent cross sections such that \t\/Q 2 < 0.1. 
No fit was made for the DD or FPD cases, while the parameter N of the profile function was 
extracted from data in the 1CDD case. 

3.2 Additional ^-dependent Do-term 

On the kinematic bin t = —0.23 GeV 2 the improvement brought by the 1CDD picture was 
significant, whereas it was less satisfactory on the bin t = —0.17 GeV 2 . It is still possible 



to add the subtraction term D of Eq. (33). This term possesses the following distinctive 
properties: 

(i) It depends on t through the multiplicative constant Tq , 

(ii) Its contribution 7i% to the CFF 7i q does not depend on £: 

ImH^fot) = 0, 



ReU q Do (U) 



n 



T(N) 



-T%{t)- 



2 1 + 2 ^AT(f + iV) 

(iii) It contributes only to the highest exponent of the polynomiality relation: 

1 1 4- (-]) n+1 T (l + -) T(N) 

ov ' s ' ' ow 2 2 ( 1+Ar ) r(2 + - + iV) 



(55) 
(56) 



(57) 



These are the properties of the D-term in the classical 2CDD formalism. 

Let us remind that the nucleon .D-term is not fixed by QCD first principles. It is however 
customary to define a flavor-singlet D-term D: 



D{a,t)= J2 ^(M> 

q=u,d,s 



(58) 



3/2. 

n 



and to project it onto the basis of Gegenbauer polynomials C-, 

oo 

D(a,t) = (l-a 2 )J2dn(t,^)C 3 n /2 ( 



a) 



(59) 



71 = 

n odd 
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Figure 5: Comparison to Hall A data : results at xb = 0.36, Q 2 = 2.3 GeV 2 (four lower plots) 
and Q 2 = 1.9 GeV 2 (upper plot), t = -0.23 GeV 2 (two lower plot) and t = -0.17 GeV 2 (three 
upper plots). The full red line corresponds to the 1CDD model, the dashed blue line to the 
classical DD Ansatz and the dash-dotted green line to the unskewed FPD limit. 



The Chiral Quark Soliton Model (xQSM) gives estimates (see Ref. |S] and references therein) 
of the first three non- vanishing terms of this expansion at a very low scale Hq ~ 600 MeV and 
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vanishing momentum transfer: 



d^ +d (t = GeV 2 ,/x 2 ) 
^ +d (£ = 0GeV 2 ,/i 2 ) 
^+ d (£ = 0GeV 2 ,/i 2 ) 



-4.0, 
-1.2, 

-0.4. 



(60) 
(61) 
(62) 



Note that Schweitzer et al. |62| report a value d\ ~ —9.46 at the low scale /Iq while Waka- 



matsu predicts d" + ~ -(4.9 - 6.2) at the same scale for the xQSM and d\ + ~ -0.716 for the 
MIT Bag model [63]. Keeping in mind the overall uncertainty on these parameters, we however 



evolve the xQSM coefficients given from Eq. (60) to Eq. (62) perturbatively at LO from the 
low scale /^ to the scale /z 2 = 2.3 GeV 2 of Hall A measurements. We follow the treatment of 
Ref. [58j EI] which assumes a vanishing gluon D-teim. at the low scale /z 2 , where both quark 
models are defined: 



rf«+ d (t = 0GeV 2 ,/i 2 ) 
+d ^-0GeV 2 ,/i 2 ) 
(t = GeV 2 ,/x 2 ) 



d u z +d (t 
dl +d 



-3.12, 

-0.71, 
-0.20. 



(63) 
(64) 
(65) 



In the specific implementation of the 1CDD formalism we are discussing here, the D-term 
is instead made of two parts. Starting from Eq. (J5l): 



D q (a) 



+l-[o| 



dPG q (P,a) 



(66) 



-i+l«l 



we restrict ourselves to the valence contribution to the DD G because the sea contribution to 
this DD is the original GK model, which is expressed in the DD representation {i.e. it is a pure 
.F-type DD). Then, from Eq. (J9l) and Eq. (32) and adding the extra ^-dependent D term (34): 



D q (a) 



a 



a 



+i-M 

i+H 
+i-M 

i+M 



d/3(/l!(/3,«) + / 9 (/3,a) 



dPfUP,<*) + D 9 (a). 



(67) 



Using the explicit expression (31) of /v al (/3,a) we find: 

«+i-M 



D q (a) 



a 



d/3 



l+|a| 



gval(/3) 




ir N ((3,a) -S(/3)7r N (0,a) 



1 7(1 -7) 1 



and define the 1CDD contribution D q CDD to the .D-term by: 

D q (a 



£) ?cdd( q; ) 



D q (a). 



+ D q (a), (68) 



(69) 
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Taking into account the support property (22) of g va i the previous equation becomes: 



D q (a) 



(i! / d/3 -^- — it N (/3, a) — air N (0,a 

o P ■> o 



a 



dp 







gval(ff) 



(1 -/5) 2 J 



-a7Tjv(0,a) / d/3 — - — 

Jl-\a\ P 



1 



2N2 2iV l j 



py 

(70) 



This D-term is indeed made of two very different contributions since only one of them depends 
on the forward limit q va \ although both contributions are of course invisible in the forward limit. 



We have seen that adding the ^-dependent Dq term (33 ) only modifies the D-term. It will not 



change the evaluation of helicity-dependent cross sections but may improve the comparison to 
helicity-independent cross sections. As the ^-dependence of Tg(£) in Eq. (70) remains unknown, 
we fit separately its values for the two different Hall A bins with \t\/Q 2 < 0.1. Since we compute 
the CFFs at LO and with GPDs evaluated at the same scale fi 2 = 2.3 GeV 2 , the fit is only 
sensitive to the total contribution D Q (t) to the charge and flavor singlet GPD H + : 



D (t) = 2(elD u (t) + eM(t) 



(71) 



where e q is the quark fractional electric charge (in units of | e |) and the sum is restricted to the 
lightest flavors because the 1CDD formalism is implemented only in the valence sector. Thus 
the parameter relevant to the fit is^j 



r (t) = 2(eX(t) + 4tf(t)). 



(72) 



At this stage it is reasonable to assume that Tq 
we fix from the data is thus: 



Tq to go further. The flavor-singlet D-term 



Diet) = 



= D u 1CDD (a) + D d 1CDD (a) + 



- D^ CDD (a) + D d 1CDD (a) + 



(r u (t) + T d (t)] 

2N2 2N 
,l + e 2 2N2^ a[ 



-a(l — a 



2\iV 



a 



2\N 

I 1 



(73) 



and is represented in Fig. [6| We compare this quantity to the yQSM estimates in Tab. [3] by 
projecting Eq. (73) onto the basis of orthogonal polynomials C n . Apart from the coefficient 
d^ +d (t = —0.23 GeV 2 ) both xQSM and 1CDD model estimates have the same order of mag- 
nitude, the latter being generically by a factor 2 smaller than the former (in absolute value). 
Remembering the small value of d\ + in the MIT Bag model reported by Wakamatsu [63j , we 
conclude that the coefficients |c^ +d | have the right order of magnitude. 
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Figure 6: D-term (73) vs x with (green dashed curve) and without (full red curve) the T 
contribution for x B = 0.36 (f ~ 0.22), £ = -0.17 GeV 2 (lower plot) and £ = -0.23 GeV 2 
(upper plot). 



7 We can also obtain this result by writing the once-subtracted dispersion relation ( 28 1 directly with the 
valence part of e^H+ + e 2 d H^ and by building the model with the corresponding PDF. 
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Coefficients 


XQSM 


Fit 


t = GeV 2 


t = -0.17 GeV 2 


t = 


-0.23 GeV 2 


dl +d 
d u b +d 


- 3.12 

- 0.71 

- 0.20 


0.39 

- 0.65 

0.12 


- 1.83 

0.018 

0.14 



Table 3: . Comparison of the coefficients of the D-term expansion Eq. (59) evaluated from the 



Chiral Quark Soliton Model and extracted from Hall A data with the Ansatz Eq. ( 70 ) 



This coefficient d\ is especially interesting because it can be related to the quark part of the 
static {i.e. defined in the Breit frame) symmetric energy momentum tensor T q 



d\{t = GeV 2 



M 



d 3 rT?Ar) 






-<fV 



(74) 



where M denotes the nucleon mass and i, j are spatial indices. Decomposing T q - in a way that 
makes manifest the distribution of pressure p{r) and shear forces s{r) of the nucleon envisioned 
as a continuous medium: 



m?) 



s q (r) 



3 



+ p q {r)5. 



i.r 



one can express d\ in terms of s g {r) or d q {r) [66J: 



d\{t = GeV 2 



M 
Y 



d fr s q {r) 



5M 



d fr p q {r). 



(75) 



(76) 



In the xQSM d™ +d {t = GeV 2 ) has a negative value, and it was conjectured from the above 
mechanical considerations that it should be the case in general. However Wakamatsu estimated 
in Ref. [63J the valence and sea contributions to the coefficient d^ + at a low scale: d\ 



l,val 



0.66 



and rf"s ea — —5.51. Both sea and valence parts may have opposite signs; this is in agreement 
with our result which applies only to the valence sector. 

In view of all the approximations involved in our comparison to the considered Hall A 
measurements (restriction to leading twist and leading order in perturbation theory, evolution 
of PDF only, etc.), we will not push further this analysis and now discuss the change in the 



fit induced by the ^-dependent Dq term (33). It allows to obtain a good agreement to Hall A 
helicity-independent cross-sections (see Fig. [Tb for the bin t = —0.17 GeV 2 . However the change 
is not dramatic concerning the bin t = —0.23 GeV 2 and the agreement is still not perfect, even 
though it is remarkable given the low numbers of fit parameters. 

From Eq. ( 69 ) we define the contribution l-i q D to the CFF 7i q associated to the D-term: 



U q {U)='H q D {U)+'H q Do {U), 



(77) 
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Figure 7: Impact of the additional £-dependent subtraction term and helicity-independent 
Hall A cross sections at x B = 0.36, Q 2 = 2.3 GeV 2 , t = -0.17 GeV 2 (above) and t = -0.23 GeV 2 
(below) for the 1CDD without .Do-term (full red line) and with Z} -t errn (dash-dotted black 
line). 



with: 



U q D {i ) t) = 2 I dxDl CDD 
Jo 



X 



U \£- x t + x 



(78) 
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and % q D (£,£) is given in Eq. (55) and Eq. (56). The contribution %£, of the D-term to the 
CFF U is thus: 



-H D = elH u D (Z,t)+e 2 d H d D (Z,t) + 



v^ T(N) T Q (t) 
2 1 + 2 ^AT(§ + iV) 2 



(79) 



since Tq (t) is extracted from a fit of D VCS data via Eq. ( 72 ) . The strength of the D-term relative 



to the the real part of the CFF 7i is given in Tab. |4j The relative weight of the subtraction term 
in ReH suggests that both components are necessary. Moreover a large D-term is required, 
i.e. a large ^-independent contribution to ReH. Such a possibility may be tested in the near 
future thanks to the forthcoming CLAS beam helicity-independent cross sections. Indeed these 
measurements will offer bins with the same t and Q 2 and allow to test the dependence of ReH, 
on£. 



£(GeV 2 ) 


ReH(T = 0) 


ReH(T ± 0) 


ReU D 


\ReU D /ReH(T o ^0)\ 


-0.17 
-0.23 


-1.59 
-1.65 


+0.20 

-1.80 


-0.52 
-2.19 


2.6 
1.2 



Table 4: Evaluation of the real part of the CFF Tt with and without the substraction term D 
and impact of the D-term (33). We remind that Q 2 = 2.3 GeV 2 and xb = 0.36. 



Even if the improvements induced by a non- vanishing T are not sufficient to get a x 2 /d.o.f. rs 
1, our study shows that changing the DD description from 2CDD to 1CDD gives results sig- 
nificantly better than what can be achieved by just adding a D-term: here the x 2 /d.o.f. goes 
down from 4.0 to 3.1 for the set of Hall A data such that \t\/Q 2 < 0.1. Meanwhile, note 



that adding a D-term such as the one in Eq. (33) to the considered 2CDD parameterization 
lowers the x 2 /d.o.f. from 5.9 to 5.1 for the same dataset. Tab. ^displays the balance between 



helicity-dependent and helicity-independent cross sections in the fits. 



Ansatz 



1CDD 
DD+D 



t = -0.17 GeV 2 



Ect 



28 / 24 
82 / 24 



Act 



158 / 48 
50 7 48 



t = -0.23 GeV 2 



Ect 



114 / 24 
445 / 24 



Act 



73 / 24 
25 / 24 



Table 5: . Statistical weight of helicity-independent (Her) and helicity-dependent (Act) cross 
sections materialized by the \ 2 P er data point for both 1CDD and DD+D Ansatze. The relative 
uncertainty on helicity-independent measurements is smaller than the corresponding helicity- 
dependent datasets, but there are many more data points for helicity-dependent cross sections. 
We considered Hall A data such that \t\/Q 2 < 0.1. 



23 



It is presumably possible to improve the agreement with the data. Indeed, the choice of 
a larger value of iV would make the 1CDD model tend to the FPD limit and thus improve 
the comparison to helicity-dependent cross sections. This would reduce the x 2 /d.o.f. because 
there are more data points for helicity-dependent than for helicity- in dependent cross sections. 
Moreover, the loss of agreement on helicity-independent cross sections may be compensated 
thanks to the additional subtraction term, but it would increase its relative weight in ReH, 
eventually putting all the GPD physics in the D-term. As said before, it will be possible to 
test this scenario on future datasets. For the time being, we consider that the other limits 
of the present phenomenological application presumably precludes from obtaining a fit with a 
X 2 /d.o.f. ~ 1. Note also the ability of the 1CDD parameterization to deal with beam helicity- 
independent cross sections, while it encounters some difficulty with beam helicity-dependent 
cross sections. Since this 1CDD parameterization was developed for spinless targets, this may 
suggest that a 1CDD picture adapted to spin-1/2 targets would improve the agreement. This 
point will be addressed in the last part of this paper. 

3.3 Comparison to CLAS data 

Beam spin asymmetries have been measured in the same kinematic range by the CLAS collab- 
oration [35J. It is natural to compare the output of the 1CDD implementation adjusted using 
Hall A cross sections to CLAS data. Considering those data in the same kinematic range, i.e. 
with t ~ —0.17 GeV 2 and \t\/Q 2 < 0.1, we compare the evaluations of beam spin asymmetries 
X 2 /d.o.f. of the 2CDD and 1CDD models. 




Figure 8: Comparison of 1CDD (full red line) and DD+D (dashed blue line) models with 
CLAS data at t ~ -0.17 GeV 2 and such that E < 0.1. From left to right: x B = 0.3205, t = 
-0.1705 GeV 2 and Q 2 = 1.9424 GeV 2 ; x B = 0.3215, t = -0.1719 GeV 2 and Q 2 = 2.217 GeV 2 ; 
x B = 0.3215, t = -0.1743 GeV 2 and Q 2 = 2.5078 GeV 2 . 



The situation is less favorable for the 1CDD implementation but can easily be explained 
by the fact that beam spin asymmetries and helicity-dependent cross sections are both mostly 
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sensitive to the imaginary part of the CFF 7i. In the 1CDD implementation IrriH is larger 
than in the FPD model, and the FPD model is in good agreement with these data. The clear 
advantage of the 1CDD parameterization appears in the description of helicity-independent 
cross sections, in particular near = 180° where helicity-dependent cross sections vanish. All 
implementations (1CDD, 2CDD and FPD) give similar estimates of helicity-dependent cross 
sections near = 60 to 90°, which is where beam spin asymmetries are maximum. The 
improvement brought by the 1CDD implementation is then hidden because there are many 
more data sensitive to IrriH than to ReH. However the flexibility of the 1CDD representation 
still exists, and we have checked that choosing N ~ 5 produces estimates of Hall A cross 
sections similar to those from the 2CDD or FPD frameworks. So even if the data do not show 
a clear advantage of the 1CDD picture, we do not loose anything by using this representation 
since it can be tuned to yield a similar quality of agreement with data. 

4 Extension to the 1CDD representation of nucleon GPDs 



4.1 From the spinless to the spin- 1/2 case 

At the time of writing this paper, Radyushkin published an extension of his earlier work on 
spinless target |32"| [H] to spin-1/2 targets [39|. The key observation is the following: it is 
natural to model H + E in the DD representation, and E in the 1CDD formalism. To see why 
it is so, let us introduce the twist-2 quark operator: 



O q 



g(0)5 7(1 iD M ...iDjg(0) 



(80) 



where D = (D — D)/2 is the covariant derivative and the operator <S projects a tensor onto its 
completely symmetric and traceless component. The matrix element of this operator between 
two nucleon states writes: 



P 



A 






O q 


2 


^ MJ,\...IMi 



p 



A 



u 



^ ' U=o ^ 



-A 



Vn-k+1 



fc=0 



+ 2_^ C n kS 



fc=0 



2M 



/•»! 



2M w 



P 



-A 



Mn-fc+1 



f^n — k 



P 



-A 



Mn-fc+1 



ftn — k 



-A, 



-A, 



-A, 



u[P 



A 
~2 



(81) 
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where the coefficients A n k, B nk and C n k depend on t. We define the DDs F q , G q and K q as the 
generating functions of these coefficients 



/?! 



(n — k) 
n ! 



(n-k)\k\J Q 



n ! 
{n-k)\k\ 



— jdf5daf5 n - k a k F q {f5 } a) = A nk 
dpda(3 n - k a k K q {(3,a) ■■ 



dpda(J n - k a k G q {f3,a) 



G n k- 



(82) 
(83) 
(84) 



The analogous of Eq. (10) in the case of a nucleon target is thus: 
A / z\ fz\ A 



P 



n-i)w K2 



2 / .. 



z 2 =0 



U\ P 



+ 



A 



7 M / d(3dae- i(i{Pz)+iai ^ 1 F q ((3,a,t) 



ia„,A u 



'fit/ 



2M 



dfida e -iP(P*)+i<x ( -¥K q {/3, a, t) 



A, 



2M 



(Az) 



d/3da e - ip{Pz)+ia * G"(/3, a, t) 



« P 



A 



+ higher twist terms. 



K) 



The GPDs P and P are defined by: 



— - e" F z ( P + — 
47T \ 2 



^~27 7 g V2 



P- 



A 

2 / z+=0 
z I =0 



H q (x,Q _f IJ A\ + / A 
2P+ l 2 / r l 2 



+P 9 (x,0« P 



A\ ia +u A, 



A 



which yields the following relations between nucleon DDs and GPDs: 



H q (x,0 



d(3da 8{x-p- af) (p 9 (/3, a) + £G 9 (/3, a) 



P 9 (x,0 = f df3da8{x-f3-aZ)(K q (f3,a)-ZG q (/3,a)\ 



It clearly appears that the DD (F + K) q generates the whole GPD (P + E) q (i.e. no P-term i 
needed in the Polyakov - Weiss gauge) and is analogous to an P-type DD in the spinless case: 
this is what we called the DD representation in Sec. 2.2 Moreover both Eq. (87) and Eq. (88) 
are reminiscent of the spinless case Eq. (J5l). We know from Sec. |1. 1| that there exists a gauge 

8 We essentially follow the presentation of Ref. [40 and use the results therein but we correct several typos. 
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86) 



(87) 
(88) 
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in which H q or E q can be cast into the 1CDD formalism. Indeed a gauge transformation of 
nucleon DDs writeqj 

F«(A«) -► ^(/3,a) + ^(/3,a), (89) 

K«{P,a) -► K*(P,a)-^(J3,a), (90) 

d\ q 
G q ((3,a) -+ G«{p,*)--gp{p,a), (91) 

and F q + K q is gauge-invariant. The extension of the 1CDD modeling of the spinless case to 
the spin- 1/2 case proposed by Radyushkin in Ref. |39| consists in modeling H + E with a single 
DD and E in the 1CDD framework. 

4.2 The regularizing role of the D-term 

Following Ref. [32l 39j , let us describe the general construction underlying the implementation 
of the 1CDD representation. In this framework GPD a GPD F generically can be expressed 
by means of a DD /: 

F(x, £) = x I d/3da 6{x-p- ag)f{P, a). (92) 

Jq 

This DD / generates a D-term: 

D(a) = a f a d 7 /(7,«). (93) 

J-i+H 

The complementary part is denoted [/]+: 

/C8,a) = [/] + G9,o)+(JC9)^-. (94) 

a 

The decomposition can be brought at the GPD level: 

F(x,0 = [F] + (x,£)+F D {x,Z), (95) 

[F]+(x,£) = x f dpda(6(x-P-a£)-6{x-atg)\f(p,a), (96) 

F D (x,0 = x [ da?^-6(x-at). (97) 

J-i « 

Let (p denote a PDF-like function (the forward limit of the GPD F) such that (f>(j3) oc l//3 a 
with < a < 1. When applying the FDDA to f(/3, a): 

f(P,a) = 7r N ((3,a)^, (98) 



3 We refer to Ref. [ID] for a rigorous treatment of the boundary conditions on DDs. 
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the D-term contribution to [F] + in Eq. (96) guarantees the convergence of the whole integral. 

smal 



Similarly to Eq. (41) and Eq. (42) we can indeed write for a valence at small x > 

-ft 



[F] + (x,0 



dp 



i-| 



- V ( ■>. — — j -7TJV (p\ - 
aA 0(/3) 







d/3 vtat /3 
ft V ?/ P 



(99) 



and: 



'" \ ( A ^— J - 7TJV f /3, 7 



iV 



4^x 






2\ JV-1 



T(2 + 2JV) 

r(i + iv) 2 



/3 + 0(/3 2 



(100) 



This is sufficient to ensure the convergence of the integral in Eq. ( 96 ) since is supposed to 



possess an integrable singularity. The general case can be derived in the same way. 

The .D-term can thus be used to regularize this 1CDD implementation and can be fixed by 
comparison to experimental or model inspired data. In such a situation this 1CDD implemen- 
tation of / is defined by the [f] + + D prescription which acts as a renormalization prescription. 

4.3 Double Distribution models for H and E 

In view of our comparison to DVCS measurements in the valence region, we still modify the 
valence sector of the GPDs H q and E q for q — u,d. We still use the GK model as a basis for 
our modifications, and describe here the modeling of the GPD E (the GPD H was described 



in Sec. 2.1). 



In the original GK model the valence part of the GPD E relies on a forward-like function 



e* al and the usual FDDA i.e. : 



Ka(*>Z,t) = [ d(3da8(x - (3 - aOM(3,uWK al (f3,t), (101) 



with: 



"val 



(/M)=p" 



-a't 



K,. 



B(l -jUval,l + ^val) 



F 



Mval 



;i-/?) 



^val 



(102) 



where B is the Euler Beta function. The values of the coefficients of Eq. (102) are given in 
Tab.! 

As announced before, we now describe H + E and E in the DD and 1CDD formalisms 
supplemented by the FDDA. Let us recall the breaking of the DD gauge invariance after the 
choice of the gauge where the Ansatz is applied. Each of the profile functions involved in either 



28 





u 


d 


K 
fJ'vaX 


1.67 

4 
0.48 


-2.03 

5.6 
0.48 


a'(GeV- 2 ) 


0.9 


0.9 



Table 6: Parameters of forward-like GPD E(x, 0, 0) in Eq. (102). 



H + E or E can come with a different parameter: 



H L ddO> + E L\ ddO> = / d/3da 6{x - p - a£)ir NlI+E (p, a)B(p) (g va i(x) + e vai (x) 






(103) 



/3 



tt Ne (j3, a) ( 5(x — /3 — a£) — 5(x — a£) 



(104) 



and the full GPD E^ is obtained after addition of the D-term D introduced to regularize the 
1CDD formalism: 



EUx,0 



HLbd(x,0+EIidb(x,0 ~ [EU+M+D 



X 



[E°J+(x,£)-D 



z 



We have chosen the following simple functional form for the .D-term: 

D(a) = Ca(l-a 2 ) } 



(105) 
(106) 

(107) 



where C depends on t and is fixed from the data, using this 1CDD implementation as a 
renormalization prescription. 

Fig. [9] shows the GPDs H™ al and _E" al with Nh+e = 1 and Ne ranging between 1 and 10. 
Since we know that the DD parameterization weakly depends on the profile function parameter, 
the model is expected to be most sensitive to the exponent of the profile function in the 1CDD 
sector. The curve corresponding to Ne = 10 is close to the FPD limit and gives an indications 
of the convergence of this family of models when Ne gets large. if" al has a characteristic 
doubly-peaked structure which is more complex than its analog in the spinless case. The peaks 
also get narrower when Ne is increased. The plot of E^ vs x oscillates less than the same plot 
of H™ al , which is understandable because H is built as the difference of two structures with an 
oscillating behavior. 
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Figure 9: Valence contributions to the GPDs H u in Eq. (105) (upper plot) and E u in Eq. (106) 
(lower plot) vs x for different values of Ne in the nucleon model for xb = 0.36, t = —0.23 GeV 2 
and Q 2 = 2.3 GeV 2 . The vertical black line signal x — £. In these plots the D-term (107) is 
arbitrarily set to 0. 



Fig. 10 compares the real and imaginary parts of the charge and flavor singlet CFF "H and E 
when varying the parameter of the profile function used in the 1CDD description for the original 
GK model, its modification inspired by the spinless case discussed in the first sections of this 
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paper, and the extension of the 1CDD model to the nucleon case. The range of possible values 
for ReH and IrnH is smaller for the spin-1/2 than for the spin-0 model. In this implementation 
of the nucleon model, the dependence of "H on N was probably softened by the addition of its 
DD part, inducing an apparent loss of flexibility. 




spin-- 1CDD 






13 



15 



N,N E 

Figure 10: Comparison of flavor and charge singlet CFFs % and S when evolving Ne (spin-1/2 
model discussed in this section) and iV (spin-0 and GK models discussed in previous sections) 
for xb = 0.36, t = —0.23 GeV 2 and Q 2 = 2.3 GeV 2 . The full red line corresponds to the spin-0 
1CDD parameterization, the dashed blue line to the DD+D parameterization, the dotted brown 
line to the spin-1/2 1CDD GPD H. 



4.4 Comparison to Jefferson Lab measurements 

We now compare all three models discussed in this paper to Hall A and CLAS data: 



DD+D GK model supplemented by a D-term as in Eq. ( 107) obtained from a fit of the data. 



spin-0 1CDD Adaptation of the 1CDD formalism in the spinless case to the GPD H as 
discussed in the beginning of the paper. 

spin-1/2 1CDD Combination of 1CDD and DD formalisms as detailed in this section. 



31 



The three models behave differently (see Fig. 11 ). The spin-0 model does a better job on beam 
helicity-independent cross sections whereas both spin-1/2 and DD+D models offer a better 
agreement with beam helicity-dependent cross sections. We note that the spin-1/2 model and 
the DD+D model are hardly distinguishable (see Tab. W\. This is also what we observe when 



comparing the three models to CLAS beam spin asymmetries (see Fig. 12). 

This surprising result may be explained by the additional D-term and the lack of sensitivity 
of these observables to GPD E. In particular, the best fit we obtained corresponds to N^+e — 1 
and Ne — > oo. From previous CFF plots Fig. [4] and Fig. 10 we see that N ~ 10 corresponds 
to the FPD limit with a good approximation, so Ne was searched between 1 and 10, and 
the fit systematically gives 10. This situation will certainly improve by adding data with a 
higher sensitivity to E, but no such DVCS measurement is available yet in the valence region. 
HERMES provide such datasets in the intermediate xb region but taking them into account 
would require the parameterization of the sea in the 1CDD framework and the treatment of a 
PDF more divergent at small nucleon momentum fractions. 



model 


t= -0.17 GeV 2 


t = -0.23 GeV 2 


DD+D 

spin-0 1CDD 

spin-1/2 1CDD 


1.9 
2.6 
1.9 


10.0 

4.1 
9.4 



Table 7: x 2 per degree of freedom for the different models and for the different values of t. 



Conclusion 

We have discussed the One and Two-Component DD representations and explained why the 
Factorized Double Distribution Ansatz breaks their equivalence. We compared these two repre- 
sentations to existing data. To keep the exercise simple, we worked with beam-helicity depen- 
dent observables in the valence region. As a first approximation, we choose to work at leading 
order, to take into account only the quark GPDs H and E, and to modify only their valence 
part. We also neglect GPD evolution and higher twist effects. Such a detailed treatment is 
beyond the scope of this paper, which simply aims at comparing the pros and cons of two DD 
representations. 

To illustrate the Two-Component DD representations (spin-0 and spin-1/2 cases) we used 
the Goloskokov - Kroll model. We built our One-Component DD models by modifying the 
valence parts of H and E and leaving the rest unchanged. Interestingly, for a given GPD, all 
models have the same unskewed limits when the profile function widths decrease to produce a 
single peak. This shows a natural limit of DD models and gives a hint about the flexibility of 
the associated DD representations. Since both One and Two-Component DD representations 
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spin-- 1CDD 
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Figure 11: Comparison to JLab Hall A helicity-dependent and independent cross sections such 
that ip < 0.1. The full red line corresponds to the spin-0 1CDD parameterization, the dashed 
blue line to the DD+D parameterization, and the dotted brown line to the spin-1/2 1CDD 
parameterization. 



have the same limit when distorting the profile function, the One-Component DD framework 
can produce similar results to the Two-component DD representation. 

The classical Two-Component DD is almost insensitive to the width of the profile function, 
while the One-Component DD displays important variations when using a spin-0 Ansatz. The 
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240 300 

<t> [deg; 



300 

<t> [deg] 



240 300 360 

4> [deg] 



Figure 12: Comparison of 1CDD (full red line) and DD+D (dashed blue line) models with 
CLAS data at t ~ -0.17 GeV 2 and such that E < 0.1. From left to right: x B = 0.3205, t = 
-0.1705 GeV 2 and Q 2 = 1.9424 GeV 2 ; x B = 0.3215, t = -0.1719 GeV 2 and Q 2 = 2.217 GeV 2 ; 
x B = 0.3215, t = -0.1743 GeV 2 and Q 2 = 2.5078 GeV 2 . The full red line corresponds to the 
spin-0 1CDD parameterization, the dashed blue line to the DD+D parameterization, and the 
dotted brown line to the spin-1/2 1CDD parameterization. 



conclusion is less clear when using a spin-1/2 modeling, which mixes H and E in the DD and 
One-Component DD representations, probably because the independent part of the model is 
not really constrained by the selected measurement sets. In that respect it would be interesting 
to carefully examine the constraints on E brought by the form factor E 2 in the spirit of the 
extensive study Ref. |67]. This point is left for future work. 

The One-Component DD has hardly been used so far in phenomenological applications 
because of its more singular behavior for small longitudinal momentum fractions. But since 
Radyushkin's treatment of divergences [32] in the One-Component DD representation, the 
implementation of the One-Component DD and Two-Component DD frameworks have the 
same complexity, at last when the GPD forward limits have integrable singularities: the same 
kind of integrals have to be dealt with. 

For these two reasons (flexibility and implementation complexity) we consider the use of 
the One-Component DD representation and its implementation along the lines described in 
Ref. [321 EH] as an interesting alternative to the earlier approach of Ref. |2"4"[ |2"5| |4"2"] . These 
features allow to build flexible ands realistic GPD models for sensitivity studies, for example to 
compute the typical size of higher-order corrections in some channels in the spirit of Ref. [58]. 

In this study GPDs were modified only in the valence region. It would be very useful 
to extend the One-Component DD implementation in a way to regularize the divergences 
associated with sea PDFs. 
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